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Abstract 

An  on-demand  single-photon  source  is  a  key  requirement  for  scaling  many  optical  quantum 
technologies.  A  promising  approach  to  realize  an  on-demand  single-photon  source  is  to  multiplex  an 
array  of  heralded  single-photon  sources  using  an  active  optical  switching  network.  However,  the 
performance  of  multiplexed  sources  is  degraded  by  photon  loss  in  the  optical  components  and  the 
non-unit  detection  efficiency  of  the  heralding  detectors.  We  provide  a  theoretical  description  of  a 
general  multiplexed  single-photon  source  with  lossy  components  and  derive  expressions  for  the 
output  probabilities  of  single-photon  emission  and  multi-photon  contamination.  We  apply  these 
expressions  to  three  specific  multiplexing  source  architectures  and  consider  their  tradeoffs  in  design 
and  performance.  To  assess  the  effect  of  lossy  components  on  near-  and  long-term  experimental  goals, 
we  simulate  the  multiplexed  sources  when  used  for  many-photon  state  generation  under  various 
amounts  of  component  loss.  We  find  that  with  a  multiplexed  source  composed  of  switches  with 
~0.2— 0.4  dB  loss  and  high  efficiency  number-resolving  detectors,  a  single-photon  source  capable  of 
efficiently  producing  20-40  photon  states  with  low  multi-photon  contamination  is  possible,  offering 
the  possibility  of  unlocking  new  classes  of  experiments  and  technologies. 


1.  Introduction 

An  on-demand  single-photon  source  is  a  key  requirement  for  many  optical  quantum  technologies,  including 
quantum  key  distribution  schemes  using  quantum  repeaters  [  1  ] ,  quantum  metrology  using  photon-number 
states  [2,3],  analog  quantum  simulators  [4] ,  and  the  boson  sampling  machine  [5] .  The  ultimate  optical 
quantum  technology  maybe  the  quantum  computer  [6] ,  capable  of  efficient  integer  factorization  and  digital 
quantum  simulation  [7, 8] ,  which  relies  critically  on  the  development  of  a  high-performance,  on- demand 
photon  source  in  order  to  efficiently  generate  large  quantum  resource  states  [9, 10] . 

While  research  efforts  continue  in  developing  on-demand  single-photon  sources  using  atomic  systems, 
achieving  high  spectral  purity  and  high  collection  efficiency  simultaneously  remains  a  challenge  [  1 1  ] .  An 
alternative  approach  commonly  used  in  quantum  optics  experiments  is  heralded  single-photon  sources 
(HSPSs)  such  as  those  based  on  spontaneous  parametric  down- conversion  (SPDC)  or  spontaneous  four-wave 
mixing:  parametric  processes  which  use  a  pump  laser  in  a  nonlinear  material  to  spontaneously  generate  photon 
pairs  (figure  1(a)).  Detecting  one  of  the  photons  using  a  single-photon  detector  will  herald  the  presence  of  the 
paired  photon  to  be  input  into  the  quantum  circuit.  Parametric  sources  are  capable  of  producing  single  photons 
with  high  spectral  purity  [11],  benefit  from  a  well-defined  wave  vector  leading  to  high  collection  efficiency  [  12] , 
and  have  no  mode-mismatch  when  integrated  on  the  same  monolithic  substrate  as  the  subsequent  circuit. 
However,  parametric  sources  are  inherently  inefficient,  due  to  the  thermal  nature  of  the  output  statistics  and  the 
requirement  to  keep  the  multi- photon  emissions  low,  resulting  in  a  maximum  single- photon  output  probability 
of  25%  [  13] .  For  applications  requiring  many  single  photons,  the  probability  of  successfully  generating  N  single 
photons  simultaneously  with  N  HSPSs  decreases  as  Nbecomes  large,  severely  limiting  the  size  of  practical 
circuits. 


©  20 1 5  IOP  Publishing  Ltd  and  Deutsche  Physikaiische  Geseiischaft 


IOP  Publishing 


New  /.  Phys.  17  (2015)  043057 


DBonneau  etal 


a) 


Pump 


O  1-4 


Pair  Source 


b)  r 

Logc  circuit 

|ES} 

;  Sn^e  photon  lQr~HI 

it 

,Tri99er  , — — — u. 

Nxl 

Figure  1.  (a)  Heralded  single-photon  source  (HSPS).  A  laser  pumps  a  parametric  pair  source  which  spontaneously  emits  a  signal  and 
an  idler  photon.  A  filter  separates  the  pump,  signal,  and  idler.  The  idler  photon  is  detected  using  a  single-photon  detector  (SPD), 
heralding  the  presence  of  the  signal  photon,  (b)  A  general  multiplexed  (MUX)  single-photon  source.  N  HSPSs  are  pumped 
simultaneously;  the  idler  photons  are  detected  while  the  signal  photons  are  stored  in  a  long  delay  line.  A  classical  logic  unit  determines 
the  configuration  for  the  N  X  1  switching  network  based  on  the  detection  signals,  routing  a  successfully  generated  single  photon  to 
the  output. 


An  approach  to  overcome  all  of  the  scaling  problems  with  HSPSs  is  the  multiplexed  (MUX)  single-photon 
source  [14] ,  which  uses  an  array  of  HSPSs,  delay  lines,  electronics  for  classical  logic  operations,  and  an  active 
optical  switching  network  to  approximate  a  true  on-demand  source  (figure  1(b)).  Using  an  array  of  HSPSs  as  a 
collective  unit  means  that  the  probability  that  at  least  one  of  the  HSPSs  emits  a  single  photon  is  high,  while  a 
switching  network  driven  by  the  heralding  signals  is  used  to  route  the  generated  photon  into  a  specific  spatial- 
temporal  output  mode.  This  results  in  the  near- deterministic  generation  of  single  photons,  allowing  for  much 
larger  quantum  circuits  than  could  be  feasibly  built  with  HSPSs  without  active  multiplexing.  MUX  sources 
inherit  the  same  benefits  as  parametric  HSPSs,  including  their  mature  theoretical  and  experimental 
investigation,  and  are  especially  appealing  due  the  prospects  of  a  fully  integrated  device  with  existing  fabrication 
processes. 

Several  theoretical  schemes  have  been  previously  investigated  [  1 5-23  ]  and  experimental  work  using  bulk 
[24]  and  integrated  [25, 26]  components  has  been  demonstrated.  Previous  work  has  highlighted  the 
fundamental  constraints  for  creating  pure  states  using  parametric  processes  assuming  ideal  MUX  components 
[  13] .  However,  in  real  physical  settings,  non-ideal  components  will  limit  the  efficiency  and  output  fidelity  of 
MUX  sources.  The  most  significant  sources  of  error  in  multiplexed  sources  are  likely  to  be  photon  loss  in  the 
optical  components  and  the  non- unit  detection  efficiency  of  the  heralding  detectors  (which  can  also  be  viewed  as 
photon  loss).  Assessing  the  suitability  of  a  MUX  sources  using  lossy  components  is  an  important  consideration 
for  many  quantum  photonic  technologies,  in  particular  the  large-scale  quantum  computer,  as  fault- tolerance 
thresholds  will  place  demands  on  the  required  source  efficiency  and  fidelity  [  10, 27] . 

The  aim  of  this  work  is  to  assess  the  impact  of  photon  loss  on  the  performance  of  multiplexed  single- photon 
sources.  We  start  by  detailing  the  heralded  single  photon  source  (HSPS),  which  is  a  key  building  block  of  all  the 
MUX  sources,  and  include  the  effect  of  photon  loss  and  non-unit  efficiency  number- resolving  and  non-number 
resolving  detectors.  We  extend  this  description  to  general  multiplexed  single- photon  sources  with  lossy 
components  and  derive  expressions  for  single-photon  and  multi-photon  emission  probabilities.  We  then  apply 
these  expressions  to  three  specific  multiplexing  source  architectures  and  consider  their  tradeoffs  in  design  and 
performance.  To  assess  the  effect  of  lossy  components  on  near-  and  long-term  experimental  goals,  we  simulate 
the  multiplexed  sources  when  used  for  many-photon  state  generation  under  various  amounts  of  component 
loss.  We  conclude  by  discussing  the  prospects  of  using  the  different  MUX  architectures  with  realistic 
components  for  near-  and  long-term  experimental  goals. 

2.  The  heralded  single  photon  source 

2.1 .  Theory  and  figures  of  merit 

All  the  MUX  sources  we  consider  are  composed  of  a  core  component  called  the  HSPS  (figure  1  (a)).  A  HSPS  is  a 
non-deterministic  source  with  a  logic  output  set  to  1  when  it  emits  a  photon  and  set  to  0  otherwise.  We  consider 
a  HSPS  composed  of  ( 1)  a  photon  pair  generation  stage  in  which  photons  are  produced  at  non-degenerate 
wavelengths  2S  (signal)  and  2;  (idler),  (2)  a  filter  which  removes  the  pump  and  separates  the  signal  from  the  idler 
into  two  different  paths,  and  (3)  a  detector  on  the  idler  arm  which  heralds  the  emission  of  the  signal  photon.  For 
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now,  we  assume  the  photon  pairs  are  generated  as  a  biphoton  two- mode  squeezed  vacuum  state  such  that  their 
joint  spectrum  is  disentangled  [28] .  Furthermore,  we  assume  that  all  signal  photons  produced  from  different 
HSPSs  are  perfectly  indistinguishable  in  all  degrees  of  freedom  except  for  the  spatial  mode  in  which  they  are 
generated. 

The  source  is  characterized  by  a  squeezing  parameter  £,  determined  by  the  pump  intensity  and  strength  of 
the  nonlinearity,  a  trigger  probability  ptrig,  the  probability  for  the  heralded  state  to  be  a  single  photon  psingle,  also 
called  the  state  fidelity  or  heralding  efficiency,  and  the  probability  for  the  heralded  state  to  be  contaminated  with 
multiple  photons,  pmuiti.  Quantifying  multi-photon  contamination  separately  from  vacuum  emissions  is 
especially  important;  vacuum  emissions  can  be  treated  as  effective  loss  in  a  linear  optical  quantum  circuit  while 
multi- photon  events  can  result  in  other  types  of  errors,  and  linear  optical  quantum  circuits  have  been  shown  to 
have  a  much  higher  tolerance  to  loss  than  to  other  errors  [  10] . 

We  aim  at  deriving  expressions  forptrig,psingie,  andpmuiti  for  the  case  in  which  the  heralding  detector  is 
number- resolving  and  the  case  in  which  the  detector  is  non-number  resolving  (also  called  a  threshold  detector). 
The  state  produced  by  the  pair  source,  assuming  spectral  disentanglement,  is  of  the  form  [29] : 


w)  =  Vi  -  m2 


|0);|0>s  +  £|l>i|l>s  + 

<  n= 2 


(i) 


where  i  and  s  are  the  idler  and  signal  modes.  In  practice,  sources  and  filters  have  losses,  and  the  heralding  detector 
does  not  have  a  unit  efficiency  detection.  We  call  rji  the  global  collection  efficiency  on  the  idler  arm  accounting 
for  all  these  effects,  and  rjs  is  the  overall  transmission  on  the  signal  arm  accounting  for  losses  in  the  sources  and 
filters.  Each  lossy  component  is  modelled  as  an  ideal  component  preceded  by  an  ideal  beamsplitter  with  a  non¬ 
unit  transmission  probability.  The  full  state,  after  accounting  for  losses  and  tracing  over  loss  modes,  can 
therefore  be  written: 


(i-m2) 


n= 0 


p=0k=0 


ii)n  -  %) 


n—k  , 


Jp,k 


(2) 


where  ppk  =  \p)i  \  k)s  ( p  \i  (k\s  and  is  the  binomial  coefficient. 

This  expression  will  be  the  starting  point  to  derive  ptrig  D,  psingle  D>  and  pmulti  D  for  the  two  detection  schemes. 
In  this  paper  we  use  the  subscript  D  on  expressions  to  indicate  the  type  of  detector  considered:  NRD  for  number¬ 
resolving  detector  and  TD  for  threshold  detector.  Detailed  derivations  can  be  found  in  appendix  A. 


2.2.  Heralded  single-photon  source  parameters 

Threshold  detector 

Starting  with  the  reduced  state  (equation  (2))  and  tracing  out  the  signal  mode,  the  probability  for  the 
detector  on  the  idler  arm  to  trigger  is  given  by  summing  all  the  contribution  of  the  states  having  at  least  one 
photon  and  is  given  by: 

Ifl 

p'rtiTD 

The  heralded  state  in  the  signal  arm  is  expressed,  renormalizing  by  dividing  by  ptrig  TD,  as: 


(3) 


P  heralded  TD 


(i  -  III2) 


Arig  TD 


Kn= 1  p=lk=0 


Tsingle  TD  ~~  (  ^  If  |2)% 


[l  -  If  I2(l  -  >7s)]2[1  -  If  l2(x  -  ^X1  -  n ;)] 


r  i*>*<n 

(4) 

) 

TD  ( 1 15  /^heralded  TD 

|l)s: 

(5) 

The  probability  that  the  heralded  state  contains  multi-photon  contamination  is  given  by 

TmultiTD  ^m=2  Is  /^heralded  TD  I  ™)sm 

_  ZTD 

P multi  TD  P single  TD  ’ 

P trig  TD 


(6) 
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with 


ZTD  =  (l  -  l£  f)  £ 
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g~\%\2  l-|l|2(l 

(i-n.) 
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(7) 


Note  that  we  can  also  similarly  derive  an  expression  for  the  probability  that  the  heralded  state  is  any  Fock  number 
state  using  (n\s  /?heralded  TD  I  n)s- 


Number-resolving  detector 

We  now  consider  the  case  in  which  the  detector  is  number  resolving.  In  practice,  we  only  need  to  discriminate 
the  single-photon  state  from  the  vacuum  and  from  states  having  two  or  more  photons.  Therefore,  the  detector 
does  not  need  to  be  able  to  distinguish  between  two  and  N  >  2  photon  states.  Starting  again  from  equation  (2), 
tracing  out  the  signal  mode  and  calculating  the  probability  to  get  only  one  photon  gives: 


Ptrig  NRD 


(i  -  III2)  III2'?,- 
( 1  —  ( 1  -*!,)  Ill2/ 


(8) 


The  corresponding  heralded  state  is: 
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As  in  the  previous  case,  we  can  compute  the  probability  the  heralded  state  contains  one  photon 

P single  NRD  ^  Is  /^heralded  NRD  |l)s: 


^single  NRD 


(l  +  (l  ->7,)(l  -ris)  III2) 

(i  -(i  -,,.)(i  -n.)  II I2)3 ; 


(10) 


Not  surprisingly,  contrary  to  the  previous  case,  number- resolving  detectors  enable  unit  fidelity  in  the  absence  of 
losses  (psingle  NRD  =  1  when  rji  =  rjs  =  1).  The  probability  for  the  heralded  state  to  contain  multi-photon 
contamination  is  given  by: 


1  -  (i -  ??i))2 

Pmulti  NRD  =  ~  "  ~  777  “  ^single  NRD  ' 


We  can  again  similarly  derive  an  expression  for  the  probability  that  the  heralded  state  is  any  Fock  number  state 
using  (n\s  /5heralded NRD  \n)s. 


2.3.  Discussion 

We  now  have  a  complete  framework  for  characterizing  any  HSPS.  We  assumed  no  joint  spectral  entanglement 
between  photons  in  a  generated  photon  pair.  However,  we  can  account  for  joint  spectral  entanglement  by 
redefining  the  fidelity  as  ps-ngle  D  =  ps ingie  D  X  P  where  P  is  the  purity  of  the  heralded  single  photon.  Recall  that, 
in  the  ideal  case  of  lossless  components,  the  probability  ptrig  NRD  of  heralding  a  unit  fidelity  single  photon  is 

bounded  by  0.25  (achieved  when  |  £  |2  =  0.5  maximizes  the  probability  |  £  |2  ^  1  —  |  £  | 2  j).  Clearly,  HSPS  do  not 
suffice  on  their  own  to  function  as  a  near- deterministic  single-photon  source,  providing  in  the  ideal  scenario  one 
photon  every  four  pulses  on  average. 

As  we  will  see  in  section  3,  the  probability  of  multi-photon  contamination  from  individual  HSPSs 
approximates  the  probability  of  multi-photon  contamination  from  MUX  sources.  We  therefore  graph  pmulti  D 
as  a  function  of  idler  transmission  rj{  for  several  squeezing  parameters  |£  |2  in  figures  2(a)  and  (b).  To  focus  on  the 
effect  of  77-,  we  take  rjs  =  1;  the  plots  therefore  serve  as  an  upper  bound  on  pmulti  D.  We  see  that  threshold 
detectors  only  reliably  herald  a  single-photon  state  for  very  low  levels  of  squeezing,  while  number- resolving 
detectors  achieve  a  much  lower  level  of  contamination  with  the  same  squeezing  parameter.  With  number¬ 
resolving  detectors,  10%  loss  in  the  idler  arm  (nz-  =  0.90),  and  pumping  with  ppair  =  0.09  ( |  £  |2  =  0.1),  the  multi¬ 
photon  contamination  level  is  almost  an  order  of  magnitude  lower  (~0.02)  compared  with  that  of  threshold 
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Figure  2.  Probability  of  multi-photon  emission  pmulti  from  a  HSPS  as  a  function  of  the  efficiency  in  the  idler  arm  rj{  for  different 
squeezing  parameters  |£  |2  using  (a)  a  threshold  detector,  (b)  a  number-resolving  detector. 


detectors  with  the  same  squeezing  parameter  (~0. 1).  Despite  their  enhanced  performance  compared  to 
threshold  detectors,  number- resolving  detectors  with  the  highest  pumping  parameters  are  not  immune  to  loss; 
even  10%  loss  in  the  idler  arm  with  ppair  =  0.25  ( |  £  |2  =  0.5)  results  in  a  probability  of  multi-photon 
contamination  of  ~0. 1 .  Achieving  a  lower  level  of  contamination  requires  a  reduction  in  either  the  squeezing 
parameter  or  the  loss  in  the  idler  arm. 

In  section  3  we  will  show  that  a  practical  multiplexed  source  will  require  operation  in  the  strong  pumping 
regime,  thus  showing  that  threshold  detectors  can  not  be  used  for  multiplexed  sources  with  the  highest  efficiency 
and  low  levels  of  multi-photon  contamination.  For  this  reason,  and  due  to  space  constraints,  we  will  only 
consider  number- resolving  detectors  in  the  remainder  of  the  paper.  However,  the  derived  expression  for 
threshold  detectors  can  still  be  used  in  the  framework  for  multiplexed  sources. 

3.  Multiplexed  single-photon  sources 

3.1.  General  considerations 

A  general  MUX  source  can  be  characterized  in  a  similar  way  to  the  HSPS.  The  probability  per  clock- cycle  that  at 
least  one  HSPS  in  an  array  of  N  HSPSs  triggers  is  given  by: 

p»r=i- <«> 

and  the  probability  per  clock-cycle  that  at  least  one  source  emits  a  triggered  single-photon  is: 

ft=w(1-(1--Ptrig)N)-  (13) 

pi  can  in  principle  be  made  arbitrarily  close  to  one,  such  that  a  near- deterministic  single-photon  source  (>99% 
emission  probability)  can  be  made  out  of  17  HSPSs  using  a  lossless  switching  network  to  route  the  photon  from 
the  HSPS  which  triggered  to  the  output  [  13] .  However,  in  any  implementation,  the  switching  network  will  have 
loss  due  to  the  optical  delay  lines — required  for  allowing  enough  time  to  reconfigure  the  switching  network 
upon  trigger  from  the  HSPS — and  the  intrinsic  loss  of  the  switch  components,  for  example  2x2  couplers  and 
phase  modulators  for  MZI-type  switches  (we  will  assume  all  2  x  2  switches  are  MZI-type  switches).  The  network 
loss,  ?7network  (iV),  is  proportional  to  the  number  of  sources  used  in  the  MUX  source,  since  additional  sources 
require  a  larger  switching  network  for  routing.  Provided  the  losses  are  equally  distributed  in  the  network  (this 
assumption  holds  for  two  of  the  three  schemes  presented  in  the  further  sections),  the  network  loss  applies  the 
same  amount  of  loss  to  every  HSPS. 

For  switching  networks  with  balanced  loss,  the  probability  per  clock- cycle  for  a  MUX  source  to  emit  a  single 
photon,  conditioned  on  the  MUX  source  triggering,  is  given  by and  is  calculated  from  equation  (5)  or  ( 10) 
(depending  on  the  type  of  detector)  byreplacing  rjs  in  the  expression  with  the  full  transmission  rjsrj  network  (N). 
Similarly,  the  probability  for  a  MUX  source  to  emit  multi-photon  contamination,  conditioned  on  the  MUX 
source  triggering,  is  given  byp^iitn  and  is  calculated  from  either  equation  (6)  or  ( 1 1 )  (depending  on  the  type  of 
detector)  by  replacing  rjs  in  the  expression  with  the  full  transmission  p5?7network  (N).  Using  this  definition  of 
Psingie,  it  then  follows  that  the  probability  per  clock-cycle  for  a  multiplexed  source  to  emit  a  triggered  single¬ 
photon  can  be  written: 
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Figure  3.  (a)  Log-tree  MUX  source.  The  trigger  output  of  each  HSPS  is  linked  to  logic  circuit  which  configures  the  setting  of  each 
switch  in  the  network,  (b),  (c)  Maximal  probability  of  triggered  single-photon  emission  and  optimal  number  of  HSPSs  per  MUX 
source  as  a  function  of  the  switch  loss.  Each  color  represents  a  given  probability  to  trigger  ptrig.  The  plain  lines  (left  axis)  represent  the 
maximal  single-photon  emission  probability  for  a  given  switch  loss.  The  dashed  line  (right  axis)  is  the  number  of  HSPSs  needed  to 
attain  the  maximal  single-photon  emission  probability  for  a  given  switch  loss,  (c)  Shows  detail  at  low  switch  losses  (<0.2  dB).  We  plot 
the  graphs  with /7delay  =  landpsingle  =  1. 


a  -  oMUX  x  nMUX 

^MUX  “  r  single  A  Ftrig  * 


For  switching  networks  with  balanced  loss,  a  convenient  lower  bound  for  this  probability  is  given  by: 

MUX  =  ^single  X  Ptrig  X  ^network  (^)> 


(14) 


(15) 


which  uses  only  the  single-photon  emission  probability  from  the  HSPSs  and  neglects  cases  in  which  multi¬ 
photon  contamination  from  the  HSPS  reduces  to  a  single- photon  due  to  loss  from  the  switching  network. 

For  a  multiplexed  source  pumped  by  a  pulsed  laser  with  repetition  rate  R ,  the  emission  rate  of  triggered 
single-photons  is: 


Ruux  -  R  X  gMUX. 

Large-scale  linear  optical  experiments  require  M multiplexed  sources  in  parallel  for  generating  M  single 
photons.  The  M-photon  state  is  heralded  by  the  simultaneous  triggering  of  all  M  MUX  sources  (at  least  one 
HSPS  per  MUX  source).  For  a  source  pumped  by  a  pulse  laser  with  repetition  rate  R ,  the  M-photon  generation 
rate  is: 


f^Mux  —  R  X  (  ^mux  )M  •  (16) 

The  multi-photon  contamination  probability  conditioned  on  heralding  the  M-photon  state  is: 

P multi  =  1  -  [l  -Pmulti  j  •  (17) 

The  analysis  so  far  has  been  kept  general,  and  can  be  applied  to  any  MUX  source  with  balanced  network  loss. 
In  the  next  sections  we  focus  on  three  specific  architectures. 

While  photon  loss  and  inefficient  detectors  are  likely  to  be  the  dominant  source  of  error  for  the  near-term 
implementation  of  MUX  sources,  we  note  that  there  are  other  sources  of  error  which  will  also  have  an  effect  on 
MUX  performance.  These  include  dark  counts  from  single-photon  detectors  [30] ,  mode-mismatch  [31],  and 
circuit  faults.  These  effects  are  left  to  a  future  analysis.  We  also  assume  that  the  detector  deadtimes  are  smaller 
than  the  system  clock  period,  and  therefore  loss  effects  due  to  detector  deadtimes  are  not  considered  in  our 
current  work. 

3.2.  Log-tree  source 

The  first  MUX  implementation  we  consider  is  called  the  log-tree  source.  In  this  scheme,  the  output  optical  ports 
of  N  HSPS  sources  are  connected  to  an  N  X  1  reconfigurable  switch,  which  is  a  logarithmic  tree  composed  of 
2x2  switches  (figure  3(a)).  Using  N sources  requires  a  log  tree  with  a  depth  of  |~  2x2  switches  in  order  to 

route  any  of  the  N  sources  to  the  output.  The  signal  photons  are  stored  in  delay  lines  as  the  electrical  trigger 
output  is  sent  to  a  logic  circuit  and  the  network  configuration  is  determined  and  set. 

The  network  loss  is  given  by  ?7network  (N)  =  rj\  ^l^deiay  where  *]  is  the  transmission  of  a  switch  in  the  log  tree 
network  and  ?/delay  is  the  loss  of  the  delay  line.  From  equation  (14),  the  probability  per  clock-cycle  for  the  log-tree 
multiplexed  source  to  emit  a  triggered  single-photon  is  given  by: 
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^tree 


and  has  the  lower  bound  (equation  (15)): 


tree 


(18) 


(19) 


The  optimal  number  of  sources  for  a  given  switch  loss  is  found  by  numerically  finding  the  N  which  maximizes 
equation  ( 19).  To  study  the  effect  of  switching  loss  in  isolation  from  other  sources  of  loss,  the  optimal  N  and  q^ee 
with  rjde]zy  =  1  and  psingle  =  1  are  plotted  as  a  function  of  the  2x2  switch  loss  in  figures  3  (b)  and  (c).  As 
expected,  we  see  that  the  probability  of  triggered  single-photon  emission  tends  towards  1  in  the  limit  of  low 
switching  loss  for  all  trigger  probabilities.  The  number  of  HSPSs  required  in  the  weak  pump  regime  (ptrig  <  0.01) 
is  likely  to  be  impractical:  obtaining  a  probability  of  single  photon  emission  q*ree  >  0.9  with  ptrig  =  0.0 1  requires 

512  sources  in  parallel  and  switches  with  0.05  dB  loss  (~0.989  transmission).  Obtaining  q*ee  >  0.9  with 
ptrig  =  0.1  requires  64  sources  in  parallel  and  switches  with  0.07  dB  loss  (~0.984  transmission).  With  the 
maximum  trigger  probability  ptrig  =  0.25,  only  16  sources  and  switches  with  0. 1  dB  loss  (~0.977  transmission) 
are  required. 

Because  the  switching  network  is  balanced,  the  probability  for  a  triggered  log-tree  MUX  source  to  emit 
multi-photon  contamination p^Sti  is  calculated  as  explained  in  section  3.1  using  ^network  (AT)  =  rj\ ^lpdeiay 
Since  the  loss  from  the  switching  network  can  only  decrease  pmuitn  the  multi-photon  contamination 
probabilities  for  the  HSPSs  in  figures  2(a)  and  (b)  serve  as  valid  upper  bounds  forp^uiti-  We  will  consider  these 
expressions  further  when  we  consider  M-photon  state  generation  in  section  4. 


3.3.  Generalized  Mach-Zehnder  interferometer 

Like  the  log- tree  source,  this  scheme  also  uses  N  HSPS  sources  connected  to  a  N  X  1  reconfigurable  switch 
(figure  4(a)).  However,  here  the  N  X  1  switch  is  a  generalized  Mach-Zehnder  interferometer  (GMZ)- 
composed  of  two  N  x  Nbalanced  splitters  enclosing  N  phase  modulators.  N  phase  modulators  are  sufficient  to 
route  any  input  to  a  given  fixed  output  port.  This  is  achieved  by  setting  half  of  the  phases  to  n  and  setting  the 
other  half  to  0,  or  by  applying  0  to  all  the  phases  to  obtain  a  full  swap.  The  NxN  passive  splitter  can  either  be  a 
N  x  NMMl  or  built  from  of  2  X  2  couplers.  Fabricating  large  N  X  Nbalanced  MMIs  with  low  loss  is  challenging, 
so  we  propose  using  cascaded  couplers  (having  a  reflectivity  of  0.5)  and  crossings  as  shown  in  figures  4(b) 
and  (c). 

The  lower  bound  on  the  probability  to  emit  a  triggered  single  photon  q^MZ  is,  from  equation  (15): 

*?GMZ  —  PsinglePdelay  —  — -Ptrig  ^  jpmodulatorPNxAP  (20) 

where  ^/moduiator  *s  transmission  of  the  modulator  section  and  rjNxN  is  the  loss  induced  by  the  balanced  NxN 

switch.  If  implemented  with  couplers  with  transmission  qcouv\er,  then  rjNxN  =  s^ow  the  effect  of 
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Figure  5.  GMZ  MUX  source:  maximal  probability  of  triggered  single-photon  emission  and  optimal  number  of  HSPSs  per  MUX 
source  as  a  function  of  the  switch  loss.  Each  color  represents  a  given  probability  to  trigger  ptrig.  The  plain  lines  (left  axis)  represent  the 
maximal  single-photon  emission  probability  for  a  given  coupler  loss.  The  dashed  line  (right  axis)  is  the  number  of  HSPSs  needed  to 
attain  the  maximal  single-photon  emission  probability  for  a  given  coupler  loss,  (c)  Shows  detail  at  low  coupler  losses  (<0.2  dB).  We 
plot  the  graphs  with  '7delay>'7modulator>  and  f>single  =  1. 


Table  1.  Comparison  of  the  requirements  and  component  depth  between  the  log  tree  and  generalized  MZI  architectures  as  a  function  of  the 
number  of  HSPSs  N.  The  first  two  columns  give  the  total  number  of  component  of  the  given  type  required  for  building  the  switch.  The  last 
two  columns  provide  the  depth  of  the  circuit  for  each  component — or  the  number  of  components  of  each  type  that  each  photon  has  to  go 
through. 


Source  type 

Total  number  of  phase  modulators 

Total  number  of  directional  couplers 

Phase  modulator  depth 

Coupler  depth 

Log  tree 

N  -  1 

2(N-  1) 

log2(N) 

2  log2(N) 

GMZ 

N 

N(jsr+iog2(Jsr)-i) 

4 

1 

2  (N  —  1) 

coupler  loss,  the  optimal  number  of  sources  to  achieve  a  probability  of  triggered  single  photon  emission  q^MZ 
W'th  Psingle'?dela>Vmodulator  =  1  is  plotted  in  figure  5. 

The  choice  between  the  log  tree  and  the  GMZ  depends  on  the  dominant  component  loss.  In  the  GMZ,  the 
photon  must  pass  only  a  single  phase  modulator,  instead  of  a  logarithmic  number  as  in  the  log- tree  scheme. 
However,  the  GMZ  requires  a  linear  scaling  in  the  number  of  directional  couplers  the  photon  passes  through, 
instead  of  the  logarithmic  scaling  given  by  the  log  tree  scheme.  Also,  the  GMZ  requires  O  (n2)  couplers  while  the 
total  number  of  components  for  the  log  tree  scales  linearly.  Table  1  shows  a  summary  of  the  different  resource 
scalings  for  the  two  architectures. 

Since  the  switching  network  is  balanced,  the  exact  probability  for  emitting  a  triggered  single-photon  gGMZ 
and  the  output  state  multi-photon  contamination  probability  pmuiti  are  calculated  as  explained  in  section  (3.1) 
usin§  ^network  =  '/delay '/modulator '/Tn'  We  wil1  consider  these  expressions  further  in  section  3. 


3.4.  Chained  sources 

In  this  scheme,  unit  cells — each  composed  of  a  HSPS,  a  delay  line,  and  a  2  X  2  switch — are  cascaded  to  build  a 
chain  of  sources  (figure  6(a)).  The  electrical  output  trigger  of  the  HSPS  directly  drives  a  switch  which  routes  the 
emitted  photon  to  the  output  port  while  the  input  port  is  routed  towards  a  blocked  path.  As  in  the  other 
schemes,  an  optical  delay  line  allows  sufficient  time  for  the  detection  of  the  idler  photon  and  the  configuration  of 
the  switch. 

The  chained  scheme  benefits  from  very  simple  control  logic  requirements,  since  each  switch  is  only  driven 
by  one  HSPS.  For  switches  driven  by  voltage,  for  example,  the  logic  consists  only  of  amplifying  the  trigger  signal 
from  the  HSPS  to  the  required  switching  voltage.  This  switching  logic  privileges  the  HSPS  emitting  closest  to  the 
output  of  the  chain,  and  hence  the  photon  from  the  HSPS  suffering  the  least  from  the  losses.  As  an  example,  if 
two  cells  are  cascaded  and  cell  2  triggers  before  cell  1,  the  switch  corresponding  to  cell  1  will  remove  the  photon 
from  the  HSPS  of  cell  2  and  replace  it  with  a  new  photon  from  the  HSPS  of  cell  1 . 

Since  the  switching  network  applies  different  amounts  of  loss  to  the  different  HSPSs,  we  cannot  derive  the 
lower  bound  for  the  single-photon  emission  probability  q*hdiin  using  the  method  given  in  section  3.1.  We 
therefore  show  a  different  derivation  here.  Given  a  chain  of  N  unit  cells,  assuming  HSPSs  with  a  probability  of 
triggering ptrig,  a  probability  for  the  heralded  state  to  be  a  single  photon  psingie,  and  switches  with  transmission  77, 
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switch  loss  (dB) 


Comparison  between  chain  (operating  at  f=1.0) 
and  tree  sources 


Figure  6.  Chain  MUX  source  and  performance,  (a)  Figure  6:  array  of  chained  sources.  N  unit  cells  are  cascaded  to  build  a  chained 
multiplexed  single-photon  source.  A  unit  cell  is  composed  of  a  HSPS,  an  optical  delay  line,  and  a  2  x  2  switch,  (b)  Number  of  cells 
required  (dash  line,  right  axis)  to  achieve  a  maximal  probability  of  triggered  single-photon  emission  q*hain  (plain  line,  left  axis)  as  a 
function  of  the  switch  loss.  We  chose  a  fractional  efficiency  of/=  0.9  to  compute  the  required  number  of  cells.  Two  regimes  are  plotted 
corresponding  to  ptrig  =  0.1  (red)  and  ptrig  =  0.25  (blue),  (c)  Comparison  between  the  chained  source  (plain)  and  log  tree  (dashed) 
triggered  single-photon  emission  probability,  as  a  function  of  switch  loss,  for  two  different  regimes,  ptrig  =  0.1  (orange)  and 
Ttrig  =  0-25  (green) .  The  optimal  single-photon  probability  q ^  is  plotted  for  both  the  chain  source — for/=  1 .0,  corresponding  to 
the  limit  of  an  infinite  number  of  cascaded  cells — and  for  the  log-tree  source — using  the  optimal  number  of  switches. 


^chain  §iven  by  summing  the  probabilities  for  each  source  to  fire  given  that  the  subsequent  ones  do  not,  and 
accounting  for  switching  losses,  which  gives: 


JS  l 

^  chain  —  ^single  ^delay  ^jP tri (  ^  _  -^trig  ) 


^single  P  delay  Ptrig  ^ 


/  \ 

~\N 

)-[ 

’] 

1  -  (!  -  -Ptrig)v 


(21) 


The  maximal  probability  for  emitting  a  single  photon  is  obtained  in  the  limit  of  an  infinite  chain 
*2max  =  Psinsie^deiay  — Practice, to  operate  at  a  fraction/of  the  maximal  probability  q^x  ,  the  length 

^Lwhich  gives: 

Ptrigj'/ 


of  the  chain  is  given  by  solving  fqimx  =  psingie1deUyptligri- 


-[( 


1  - 


■(- 


N  = 


In  ( 1  —  x) 

ln[(1-rtrig)'/] 


(22) 


We  show  the  single-photon  emission  probability  q*hain  with  psingle  ?7delay  =  1  and  the  number  of  cells 
required  N  as  a  function  of  the  switching  loss  for/=  0.9  in  figure  6(b).  As  with  the  log  tree  architecture,  the 
required  number  of  sources  becomes  impractical  for  small  ptrig,  but  can  be  kept  below  8  when  operating  at 
ptrig  =  0.25.  The  performance  of  the  chained  scheme  compares  well  with  the  performance  of  the  log  tree,  as 
shown  in  figure  6(c).  It  has  better  performances  for  ptrig  =  0.1  (for  a  switch  loss  less  than  1.6  dB)  and  is  close  to 
the  log- tree  performances  for  ptrig  =  0.25. 

The  derivation  for  the  exact  probability  of  triggered  single  photon  emission  qc hain  and  the  output  state  multi¬ 
photon  contamination  probability pmuiS  are  shown  in  appendix  B. 
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a)  b) 


Figure  7.  Left  axis  (plain  lines):  requirements  for  generating M  single  photons  at  a  rate  of  100  Hz,  from  a  100  MHz  pulsed  laser  for  two 
different  idler  detection  efficiencies.  The  red  lines  correspond  to  a  lumped  detection  efficiency  of  rji  =  rjs  =  0.9.  The  blue  lines 
correspond  to  a  lumped  detection  efficiency  of  rji  =  rjs  =  0.99.  (a)  Maximum  tolerable  switch  loss  (per  2  x  2  switch)  versus  number 
M  of  log-tree  MUX  sources  operated  in  parallel.  Right  axis  (dashed  lines):  probability  of  multi-photon  contamination  in  output  state 
from  the  M  MUX  sources  when  using  the  maximum  tolerable  switch  loss,  (b)  Number  of  HSPSs  required  per  MUX  source  in  order  to 
meet  the  requirements  in  (a). 


4.  Discussion:  generating  M  single  photons  using  M  multiplexed  sources 

Experiments  with  about  eight  single  photons  represent  the  current  state-of-the-art  photon  number  for  SPDC 
experiments  in  practice,  so  we  will  examine  the  performance  of  multiplexed  sources  for  the  generation 
of  >10  single  photons  in  separate  spatial  modes.  We  first  consider  using  M log- tree  multiplexed  sources  to 
generate  M  single  photons  at  a  rate  of  100  Hz,  assuming  a  pulsed  laser  seed  with  a  repetition  rate  of  100  MHz, 
number-resolving  detectors,  and  ppair  =  0.1.  Figure  7  shows  the  maximum  tolerable  switch  loss  for  generating 
M  single  photons  under  these  assumptions,  and  the  resulting  probability  of  multi-photon  contamination  in  the 
output  state  for  this  same  switch  loss.  We  considered  two  values  for  the  lumped  idler  detection  efficiency  (which 
includes  filter  loss),  rji  =  0.9  and  rji  =  0.99,  and  assumed  that  this  loss  is  also  applied  to  the  signal  channel.  We 
see  that  with  rji  =  0.9  and  switches  with  1  dB  loss  (~0.794  transmission),  14  single  photons  can  be  generated  at 
the  target  100  Hz  rate  with  multi-photon  contamination  below  10%  .  Extending  beyond  this  regime  to  perform 
experiments  with  20-40  photons  with  less  than  10%  multi-photon  contamination  requires  rji  «  0.99  lumped 
idler  detection  efficiency  with  2  x  2  switches  having  each  <0.4  —  0.2  dB  loss  (~0.912  —  0.955  transmission). 

Although  these  numbers  are  beyond  what  is  attainable  using  the  current  state-of-the-art  components,  they 
maybe  attainable  in  the  near-to-mid-future  with  further  development  and  integration.  For  reference,  prototype 
all-optical  switches  in  fiber  have  demonstrated  insertion  loss  of  ~0.6  dB  (~0.871  transmission)  [32] ,  on-chip 
filters  for  sufficient  pump  suppression  have  demonstrated  insertion  loss  of  of  ~0.8  dB  (~0.832  transmission) 
[33] ,  and  state-of-the-art  single-photon  detectors  have  demonstrated  efficiencies  of  0.93  [34] .  Although  the 
single-photon  detectors  referenced  in  [34]  are  non -number- resolving  and  are  limited  to  a  deadtime  of  40  ns,  we 
note  that  passive  multiplexing  techniques  using  arrays  of  detectors  [35, 36]  can  achieve  approximate  number¬ 
resolving  capabilities  and  reduce  deadtime.  However,  the  exact  effect  of  loss  using  these  specific  detector 
architectures  on  multiplexed  sources,  and  their  resource  requirements,  is  beyond  the  scope  of  this  paper. 

For  applications  requiring  a  very  large  number  of  single-photons  (M  >>  40),  such  as  the  universal  quantum 
computer,  error- correcting  and  loss-tolerant  encodings  can  be  used  to  achieve  fault- tolerant  computation, 
provided  the  error  rates  are  below  required  thresholds.  Meeting  these  thresholds  will  require  sources  with  the 
highest  single-photon  emission  probability  and  very  low  multi-photon  contamination.  We  briefly  consider 
what  is  attainable  with  optimistic  long-term  goals  for  the  performance  of  required  hardware  components. 
Assuming  0.99  lumped  idler  detection  efficiency,  2x2  switches  with  0.98  transmission,  and  ppair  =  0.1,  this 
corresponds  to  a  multiplexed  single-photon  source  efficiency  of 0.8744,  a  multi-photon  emission  probability  of 
0.0017,  and  64  HSPSs  per  MUX  source.  Using  ppair  =  0.25,  this  increases  the  single-photon  source  efficiency  to 
0.8965  and  reduces  the  resource  requirements  to  16  HSPS  per  MUX  source,  although  at  the  cost  of  an  increased 
multi-photon  contamination  probability  of 0.0083.  Further  research  is  necessary  to  translate  these  single-  and 
multi-photon  emission  probabilities  into  error  rates  for  fault- tolerance  in  linear  optical  quantum  computing, 
and  to  optimize  multiplexed  source  architectures  for  scalable  quantum  computation. 
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5.  Conclusion 

In  our  analysis  we  have  studied  multiplexed  single-photon  sources  with  photon  loss  and  inefficient  detectors. 

We  derived  expressions  for  single-photon  emission  probability  and  multi-photon  contamination  probability 
using  number- resolving  and  non- number- resolving  detectors  for  a  general  HSPS  and  for  three  MUX 
architectures:  a  log- tree  switching  scheme,  a  chained  switching  scheme,  and  a  generalized  MZI  scheme. 

Our  findings  indicate  that  number- resolving  detectors  offer  a  considerable  advantage  over  threshold 
detectors  for  MUX  sources,  and  are  essential  for  MUX  sources  with  the  highest  efficiency.  The  performance  of 
the  chained  scheme  compares  well  with  the  performance  of  the  log  tree,  and  has  a  simple  logic  requirements  that 
maybe  beneficial  for  a  near-term  implementation.  The  GMZ  offers  an  alternative  switching  structure  that 
requires  fewer  phase  modulators  but  more  directional  couplers  than  the  other  schemes. 

All  three  of  the  architectures  considered  are  capable  of  a  high  efficiency  single-photon  source  with  low  multi¬ 
photon  contamination  if  the  components  are  close  to  ideal.  For  example,  with  a  log- tree  architecture  composed 
of  switches  with  <0.4  -  0.2  dB  loss  (~0.912— 0.955  transmission)  and  number- resolving  detectors  with  99% 
efficiency,  a  single-photon  source  capable  of  producing  a  100  Hz  rate  of 20-40  photon  states  with  less  than  10% 
multi-photon  contamination  is  possible.  Such  a  source  would  be  a  valuable  resource  for  quantum  technologies 
such  as  the  boson  sampling  machine  and  quantum  simulators.  With  lower-loss  switches  and  the  high  efficiency 
detectors,  a  MUX  source  approaching  the  threshold  requirements  for  a  fully  fault-tolerant  universal  quantum 
computer  should  be  possible.  Further  work  is  necessary  to  find  the  optimal  multiplexing  schemes  which  could 
encompass  the  direct  generation  and  multiplexing  of  multi-photon  states  or  entangled  resource  states. 
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Appendix  A.  Heralded  single-photon  source  parameters  derivation 


Threshold  detector 

The  source  output  state  after  accounting  for  losses  is  given  by: 


oo  n  n 


\n—k  . 


(23) 


p = t1  -  i£f)  Pcw('  - u  pp* 

v  n=0  p=ok=0 

where  pp  k  =  \p)j  \  k)s  ( p  |z-  (k  \s  and  Cj)  is  the  binomial  coefficient. 

The  probability  for  the  detector  placed  on  the  idler  arm  to  trigger  is  given  by  summing  all  the  contribution  of 
states  having  at  least  one  photon.  Starting  with  the  reduced  state  after  tracing  out  the  signal  mode: 
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We  can  check  the  validity  of  the  expression  in  the  two  extreme  cases.  If  r\{  =  1,  the  probability  to  trigger  is 
maximal,  plllg  =  |£|2.Ifp,  =  0  then  ptrig  =  0. 

The  heralded  state  in  the  signal  arm  is  expressed,  renormalizing  by  dividing  by  ptrig,  as: 


oo  n  n 


i>„  n, = f>  I,  2 1«  f  2  2 w ( 1  -  ( 1  -  i.r V 

T trig  TI)  m=  1  \  n=0  p=0k=0 
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Ttrig  TD 


oo  n  n 
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v  n=l  p=lk=0 


We  can  compute  the  probability  that  the  heralded  state  is  a  single  photon  using  psingleTD  =  (l|s  pheralded  XD  |l)s: 
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Using  the  identity  (obtained  from  taking  the  derivative  of  the  geometric  series)  =  l/^  —  x)2; 

i  t1-^) 
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[i  -  II I2 ( i  -  n s)]  [l  -  |||2(i  -  ps)(i  -  //,•)] 

In  the  lossless  case,  we  obtain  psingle  TD  =  (l  —  |^|2j,  which,  with  a  triggering  probability  ptrig  =  |£|2(from 

equation  (24)),  is  consistent  with  a  probability  of  getting  one  photon  of  |  £  |2  ^  1  —  |£  |2)  (from  equation  (1)). 

The  probability  that  the  heralded  state  contains  multi-photon  contamination  is  derived  as  follows.  First,  we 
define: 
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Then  it  follows  that  PmultiTD  =  £”=2<m|s  ^heralded  TD  Ms- 
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Number-resolving  detector 

Now  we  consider  the  number- resolving  detector.  Starting  with  equation  (2),  tracing  out  the  signal  mode  and 
calculating  the  probability  to  get  only  one  photon: 
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We  can  check  the  validity  of  the  expression  in  the  two  extreme  cases.  If  rji  =  1,  the  probability  to  trigger  is 
maximal,  ptrig  =  (l  -  |||2)  |£|2.  If  ^=0  then  ptrig  =  0. 

The  corresponding  heralded  state  is: 
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As  in  the  previous  case,  We  can  compute  the  probability  that  the  heralded  state  is  a  single  photon  using 
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The  probability  that  the  heralded  state  contains  multi-photon  contamination  is  derived  as  follows.  First,  we  define: 
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Then  it  follows  that  the  probability  that  the  heralded  state  contains  multi-photon  contamination  is  given  by 

TmultiNRD  =  Hm=2^m  Is  /^heralded  NRD  \m)s: 
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Appendix  B.  Chain  source  derivation 

The  exact  probability  of  single-photon  emission  qc  hain  can  be  calculated  as  follows.  The  state  from  a  chain  of  N 
sources  is: 


N—l 

Pd  ~  -Ptrig  D  ^  ^  -PtrigD^  PjD » 

j=  0 


where  pj  D  is  the  state  heralded  by  an  individual  HSPS  for  which  the  signal  arm  goes  through  j  +  1  switches  each 
having  transmission  rj ,  and  D  is  the  type  of  detector  (TD  for  threshold  detector  or  NRD  for  number  resolving 
detector). 

For  a  number  resolving  detector: 
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The  probability  that  the  MUX  source  triggers  is  given  by: 

The  probability  of  single-photon  emission  conditioned  on  the  MUX  source  triggering  is  given  by: 

Ps^D  =  -4^(1\sPd  |1>. 

”  trig 

1  N_1  j 

~~  chain  2^trig  D  (  ^  —  P trig  D  )  P single;  D > 
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where  psingle D  is  the  probability  of  single-photon  emission  for  an  individual  HSPS  heralding  the  state  p ■  D 
(calculated  in  the  same  way  as  in  section  2  but  using  pjD).  Then  we  have: 


_  ^  chain  v  ^  chain 

^  chain  D  -^single  D  PtrigD’ 


The  probability  of  multi  photon  contamination  conditioned  on  the  MUX  source  triggering  is  given  by: 


P2m  = 
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with 


Pmulti;D=  P j  D  \m)s 

m= 2 

=  Zj  £>  -  psingle  .  D , 

where  Z;  D  is  the  normalization  constant  for  an  individual  HSPS  heralding  the  state  pj  D  (calculated  in  the  same 
way  as  in  section  2  but  using  pj  D)  Then  we  can  rewrite  (29)  as: 
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where 
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